Unconventional superconductivity on honeycomb lattice: the theory of Kekule order 
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A spatially non-uniform superconducting phase is proposed as the electronic variational ground 
state for the attractive interactions between nearest neighbors on graphene's honeycomb lattice, close 
to and right at the filling one half. The state spontaneously breaks the translational invariance of the 
lattice into the Kekule pattern of bond order parameters, and it is gapped, spin triplet, and odd under 
the sublattice exchange. With the increase of attractive interactions we first find the transition from 
the semimetallic phase into the p-Kekule superconductor, defined as being odd under the exchange 
of Dirac points, with the additional discontinuous superconductor-superconductor transition into 
the even s-Kekule state, deep within the superconducting phase. Topological excitations of the 
Kekule superconductor and its competition with other superconducting states on the honeycomb 
lattice are discussed. 
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I. INTRODUCTION 

Fermions on graphene's honeycomb lattice can in prin- 
ciple find themselves in a plethora of insulating phases, 
depending on the relative magnitudes of different compo- 
nents of a finite-range repulsive interaction, for example 
[l| . If the net interaction would have an attractive com- 
ponent, on the other hand, there would be a variety of 
superconducting states available to Dirac quasiparticles 
for pairing and condensation. Some of them are quite 
conventional: the on-site attraction would clearly favor 
the usual s-wave singlet pairing [2|. Others are already 
less so; the second-nearest-neighbor attraction, for exam- 
ple, leads to an f-wave superconductor Q , which changes 
sign six times around the Brillouin zone. Another exotic 
superconducting state on honeycomb lattice was argued 
to arise from the nearest- neighbor attraction Q : instead 
of gapping the Dirac points it lowers the energy of the 
Dirac-Fermi sea by effectively increasing the Fermi ve- 
locity. Only away from half-filling does this state acquire 
a finite superconducting gap, which is then proportional 
to the chemical potential. A closely related supercon- 
ducting ground state was also discussed in the context of 
the t-J-U model and graphite [5[- This hidden supercon- 
ducting order is otherwise a spin-singlet, and even un- 
der the exchange of the two sublattices and/or the Dirac 
points. Since the electrons in graphene have three sets 
of discrete indices, the sublattice, valley, and real spin, 
possible superconducting states may exhibit various sym- 
metries with respect to spatial and time inversions [6|, |7| . 
Together with the observation of superconductivity in 
graphite [3, the intricate structure of the superconduct- 
ig vortex [Tl |9|, llO| , novel proximity phenomena ll| and 



the quantum criticality |12l Il3l | , this makes the problem 
of superconductivity in graphene or in an optical honey- 
comb lattice engaging from theoretical as well as experi- 
mental points of view. 

In this paper we will be concerned with the forms of the 
superconducting condensate on the honeycomb lattice 
at, and therefore also near, half-filling. As a convenient 
point of departure we consider the problem of graphene 



with the chemical potential right at the Dirac point and 
with the attraction only between the electrons residing 
on the nearest neighbors of the honeycomb lattice. The 
motivation for studying such a pairing interaction of a 
finite range comes in part from the theories of boson- 
fermion mixtures in optical lattices, where the nearest- 
neighbor attraction between fermions arises upon inte- 
gration over the bosonic degrees of freedom [IJ]. Also, 
since the fermions in reality certainly experience a strong 
repulsion when they find themselves on the same site, the 
attraction between the nearest neighbors appears to be 
the simplest reasonable assumption that would still lead 
to pairing. Our conclusion about the superconducting 
ground state that arises as the ECS mean-field solution in 
this model is unusual and qualitatively different from the 
previous study [43. Within the standard mean- field ap- 
proach we find the superconducting state with the lowest 
energy to be the spin-triplet, non-uniform condensate, 
which is odd under the exchange of the two sublattices. 
The spatial Kekule pattern [15| of bonds between the 
paired electrons on nearest-neighbors has the periodicity 
of 2(5, where ibQ are the Dirac points, which allows it 
to connect the two Dirac valleys and that way open the 
mass-gap in the Bogoliubov quasiparticle spectrum. It 
is an example of Fulde-Ferrell-Larkin-Ovchinikov [la. [l7| 
type of superconducting phase appropriate to the honey- 
comb lattice. 

We argue that the development of such a non-uniform 
superconductor at T=0 may preempt the formation of 
the previously proposed hidden order, which in our ap- 
proximation we indeed find to be suppressed at all cou- 
plings. The Kekule superconductor breaks the exact 
particle-number and the spin-rotational symmetries, and 
exhibits three massless and three massive modes in the 
ordered phase. It also rather weakly breaks the internal 
and approximate U{\) symmetry between various Kekule 
patterns. Near the semimetal-superconductor transition 
we find the p-Kekule state, odd under the valley ex- 
change, to have the lowest energy, with an additional dis- 
continuous transition within the superconducting phase 
into the s-Kekule state, even under the valley exchange. 




FIG. 1: The unit cell of the Kekule lattice of superconducting 
bond order parameters. The red line corresponds to A cos a, 
the bold line to Acos(a+27r/3), and the thin line to Acos(a — 
27r/3). The unit cell contains six sites (blue points) and nine 
bonds. When periodically arranged in a triangular lattice of 
period 3 it yields the Kekule pattern. 



at a stronger attractive interaction. 

The target space of the Kekule order parameter is 5*3, 
the surface of sphere in four dimensions. The topology 
of this space implies that there are no stable topologi- 
cal defects in our two-dimensional system, and therefore 
presumably no sharp finite temperature phase transition. 
Explicit breaking of the rotational symmetry, by an ex- 
ternal magnetic field or the spin-orbit interaction, for in- 
stance, changes the target space for the order parameter 
and restores the possibility of topologically distinct de- 
fects. The cases of easy plane and easy axis, introduced 
by the two terms mentioned above, are both discussed. 
We also list all other gapped and hidden (gapless) super- 
conducting states on the honeycomb lattice, and briefly 
comment on their competition. 

The paper is organized as follows. In the next sec- 
tion we write the Bogoliubov-de Gennes Hamiltonian in 
the Dirac form, and introduce the non-uniform Kekule 
ansatz for the superconducting bond order parameters. 
The minimization of the mean-field energy for the sim- 
plest non-uniform state and the resulting s-Kekule su- 
perconductor is presented in section III. In sec. IV we 
discuss the competition between the Kekule and hidden 
orders. In sec. V a more general Kekule pattern is con- 
sidered, and the p-Kekule state is defined. The other 
possible superconducting orders are discussed in sec. VI, 
and the issue of topological defects and the target space 
for the Kekule order parameter in sec. VII. Concluding 
remarks are given in sec. VIII. 



II. BDG-DIRAC HAMILTONIAN AND THE 
KEKULE ANSATZ 



Consider the usual tight-binding Hamiltonian for spin- 
1/2 fermions on honeycomb lattice at half-filling, with an 



attractive interaction between the nearest neighbors, 

H = Ht-vY, E M^W'iy), (1) 

{x,yj cr,(T'=t4 



Ft = <^ ^ u^(xK(y) + /i.c., (2) 

{x,y) cr=t4 

where V > 0. Uaix) and Va{y) are the fermionic opera- 
tors at the two triangular sublattices of the honeycomb 
lattice. Decoupling the interaction term in the particle- 
particle channel yields the Bogoliubov-de Gennes (BdG) 
Hamiltonian 

HsdG^ Ht-^ E l^aa'{y,x)ul,{x)vl{y) + h.c., 

{x,y} a,a'=t,l 

(3) 
with the superconducting order parameters to be deter- 
mined self-consistently as 



^aa'iy^x) = V{v„{y)u„'{x)). 



(4) 



Assuming the order parameters to be much smaller 
than the bandwidth the condensation energy comes 
mainly from the pairing of the quasiparticle states near 
the two Dirac points. Let us form a 16-component Dirac- 
Nambu fermion * = (*p, '^hV , with *p = (*p^, *p;)^ 
and *,, = {^hu -*/it)^, and 

^IM = {ua{Q + q),v„{Q + q),u,{-Q + q),v,{-Q + q)), 

(5) 



*L('?) - {vl{Q~q),ul{Q-q),vl{-Q-q),vl{-Q-q)). 

(6) 
The tight-binding Hamiltonian at low energies then be- 
comes 



Ht=Y,^\q)HD^{q) + 0{q^), 



(7) 



with Hd as the Dirac Hamiltonian in two dimensions, 
which in our construction and in the first quantization 
assumes a particularly simple form. 



HD = TQ®aQ® ijojtqi 



(8) 



Here, 70 = ctq 0-3, 71 = 0-3 (^ a2, 72 = ctq 18) cti , are the 
usual four-component anticommuting Hermitian gamma- 
matrices [1[. The two-component Pauli matrices {tqjt} 
operate on Nambu's, and {cto,!?} on the spin indices. 
We will also define the remaining two gamma-matrices 
as 73 = CTi (81 172, and 75 = (T2 (8) cr2- For convenience, 
hereafter we also set the Fermi velocity vp = \/3t/2 — 1 
and the lattice spacing a to unity. 

Next, we define the Kekule ansatz for the supercon- 
ducting order parameter: 

A„„{x,y) ^ A„cos{Q- {x + y)+a), (9) 



-{A:^ti^,y) + AuiS,y}) ^ AcosiQ ■ (x + y) + a), (10) 



the above ansatz the BdG Hamihonian can be rewritten 
as 



-iA^^{x,y)-A^^{x,y})^A'. 



(11) 



The components of the triplet are assumed to be spatiahy 
periodic, with the periodicity of 2Q, whereas the singlet 
component is simply uniform. The "angle" a parame- 
terizes different spatial patterns of the order parameter. 
The unit cell of the Kekule lattice is depicted in Fig. 1. 



HbcIG = Ht+J2 ^^©(M + Af )*('Z) (12) 



where the two matrices appearing in the last term are 



III. S-KEKULE GROUND STATE 



M' = i{Re[A']Ti + Im[A']T2) ® ctq ® ilol?,)HD, (13) 



We determine first the optimal Kekule ground state for 
a — 0, and then consider a more general solution. With 



and 



J 



M = [{R+T2+I+Tl)®(J2 



{I-T2 - R-Ti) (8) (71 + {Xti - Yt2) (g) 0-3] (Ki 70, 

I 



(14) 



where 



R± 



-(i?e(At)±i?e(A^)), 



and 



/± = -(/m(A^)±/m(AJ), 



A = X + iY. 



(15) 



(16) 



(17) 



Before proceeding with the diagonalization of the BdG 
Hamiltonian it is worth pausing to register its sym- 
metries. The Dirac Hamiltonian commutes with N = 
^3 (8) ctq <8) / and P = T3 (g) ctq <8) 17375 , which in our rep- 
resentation stand for the particle-number operator and 
the generator of translations. It also commutes with 
Ik = T'o'X'croig) 17175, and /„„ = tq (g) ctq ^ 72 , when accom- 
panied with the axis inversions qi -^ —Qi, and q2 -^ —92, 
respectively. The latter two operations represent the ex- 
changes of the two Dirac-points and the two sublattices, 
respectively [l|. Hd also commutes with all three gen- 
erators of rotations of electron spin, S = tq (E) (J ® I- 
The matrix M does not commute with iV, P, and /„„, 
but, for a ~ under consideration at the moment, it 
does commute with Ik- It therefore represents a spa- 
tially non-uniform superconducting condensate, which is 
odd under the sublattice exchange and even under the 
exchange of Dirac points. We will call it the s-Kekule su- 
perconductor. Since it violates the spin-rotational sym- 
metry, the matrix Af represents a triplet superconducting 
state, which breaks two generators of spin rotations. 



The matrix M' , on the other hand, is a product of the 
Dirac Hamiltonian and another matrix which, in our rep- 
resentation, by itself would represent the singlet s-wave 
order parameter jl8{ . Since the two factors anticommute 
the presence of the imaginary unit in Eq. (13) makes 
the matrix M' Hermitian. The matrix M' represents the 
hidden superconducting order [4|. This superconducting 
state, however, suffers from an energetic disadvantage: 
since A/' vanishes precisely at the Dirac points, opening 
of the order parameter A' seems like an ineffective way 
to lower the energy of the filled Dirac-Fermi sea. We 
will argue shortly that in competition with the Kekule 
triplet the hidden order is likely to be energetically in- 
ferior. Therefore we set A' = for the time being, to 
return to the issue of the hidden superconducting order 
only after we determine the optimal triplet s-Kekule or- 
der parameters. 

Setting then A' = one finds 

[Ho+Mf = {q^+m^){To(g)ao(g>I)+2T3(g){n-a)(g)I, (18) 
where the vector ft has the components: 

ri = {XR+ + YI+,YR_ - XI.,R+R_ + /+/_), (19) 

and the mass-gap is 

m^ =X'^ + Y'^ +R\+Il+R^_+I'i. (20) 

The mean-field ground state energy per site of honey- 
comb lattice is therefore 



E _ 3m2 ^ 



s=± 



dq 



iq^+m^ + 2s\n\y/^, (21) 



where N is the number of points in the first Brillouin 
zone. We also assume an ultraviolet cutoff A in the in- 
tegral over momenta, which is here performed only near 
the two Dirac points. Differentiating with respect to \n\ 
immediately shows that for any value of the mass m the 
minimum of energy lies at \n\ = 0. We set therefore 
ni = and 712 = 0. Viewed as a set of two linear equa- 
tions for the variables X and Y they will have a non- 
trivial solution only if 



R^R- 



Ul- = 0, 



(22) 



which also happens to be the remaining equation n^ — 0. 
The trivial solution X = Y = will be discussed sepa- 
rately in sec. VII. The condition |n| = yields therefore 
only two, and not three independent equations. The last 
equation then implies 



I A. 



tl 



HI 



(23) 



and the remaining condition constrains the order param- 
eter's phases as 



<p^ + (pl = 2(j) + TT, 



(24) 



where Ag- = \Aa-\exp{i(pcr), and A = |A| exp(i(/)). Fi- 
nally, the minimum of energy is at the value of the mass- 
gap m = Too determined by the gap equation: 






2V 
3 



dq 



(27r)2 {q^+miy/^' 



(25) 



which has a solution for V > Vc- At the minimum, af- 
ter some straightforward algebra the matrix M can be 
written as: 



M — too(ti COS0 — T2 sine/)) (g) [sin0((Ti cos((/)4, 

r 



-I- 0-2 sin((/); -(/))) -I- 0-3 COS 6*] (g) 70, 



(26) 



where |A| = mo cos 9, and |A^| — |A^| = TO,osin6'. At the 
minimum of the energy the Kekule state has three hard 
and three soft modes: the angles {9, (/>4. — 0) determine the 
preferred spin axis for the triplet state, and is the super- 
conducting phase. Note that the condition for the energy 
minimum |n| = eliminated three out of six linearly in- 
dependent matrices that appear in the matrix M in Eq. 
(14). The remaining three matrices anticommute among 
themselves as well as with the Dirac Hamiltonian and 
therefore enter as a sum of squares into the expression of 
the ground state energy. This quite generally appears to 
be the optimal way for the filled Dirac- Fermi sea to lower 
its energy. Another example of this rule is the emergence 
of the easy plane for the Neel order parameter for the 
antiferromagnetic state on the honeycomb lattice in the 
magnetic field [19]. Further consequences of this rule for 
the form of the order parameter in the presence of the 
terms that break rotational symmetry will be discussed 
in sec. VII. 



IV. HIDDEN ORDER PARAMETER 

Let us now restore the possibility of the hidden super- 
conducting order, while retaining the energy-minimum 
condition n = 0. The mean-field energy per site is now 
modified into: 



E 

2N 



3(to2 



2|A' 



2V 



dq 



(27) 



In writing this expression we assumed the relative phase 
between the hidden and the Kekule order parameters to 
be 7r/2, so that the matrices M' in Eq. (13) and M 
in Eq. (26) anticommute and enter the energy expres- 
sion as a sum of squares. The relative factor of two in 
the first term derives from the sum of order parameters 
over a Kekule unit cell (Fig. 1). The critical interac- 
tion for the appearance of the s-Kekule order is therefore 
Vc — 37r/A, whereas for the hidden order, in absence of 
the Kekule state, it would be V^ = IStt/A^. Choosing the 
cutoff A even as big as unity, which, for instance, would 
represent the interval of the energies over which the tight- 
binding density of states is approximately linear, we see 
that by increasing the interaction at V = Vc the system 
first becomes the Kekule superconductor, with ttiq ^ 0. 
Upon further increase of the interaction the amplitude of 
the order parameter to,o grows, and then suppresses any 
appearance of the hidden order. We believe the reason 
for this outcome of the competition to be quite physical: 
given the choice whether to open the gap in spectrum 
or increase the velocity of excitations, all the rest be- 
ing equal, the system chooses the former option as ener- 
getically preferable. The reader should be warned, how- 
ever, that this conclusion could in principle be overturned 
upon inclusion of the states farther from the Fermi level 
into the energy calculation. The pure hidden order, or 
even the coexistence of the two orders, seem conceivable 
as well. Since the presence of the residual repulsive inter- 
actions in a real system will always broaden the single- 
particle states away from the Fermi level, it is difficult 
to say anything more definite on this issue beyond the 



low-energy approximation we employed. 

We have checked, nevertheless, that our conclusion re- 
mains unaltered within the present mean-field calcula- 
tion that keeps all quasiparticle states perfectly sharp, 
upon inclusion of the states from the entire first Bril- 
louin zone. This way we find the two critical interac- 
tions defined above to be VJ = 3/0.786 = 3.816 and 
Vc = ((3/2)/0.727) = 2.063, in qualitative agreement 
with the conclusion based on the linear approximation 
to quasiparticle dispersion. For further details of this 
computation the reader should consult the Appendix. 




V. P-KEKULE STATE 

We turn to a general Kekule state with the parameter 
a y^ next. Select the spin axis so that A-j- = A^ = 
by setting the angle 9 = in Eq. (26). Without a loss in 
generality one may choose then the order parameter A 
to be real, and write the BdG Hamiltonian in real space 
as 



H 



BdG — 7 , 



<^Hxi)[{Ta (g,T)+ A(ti (g, K)]^{x2), (28) 



where xi and X2 belong to the same sublattice, and 

$^(x) = {u^{x),vi^{x + b),ul{x),-vl{x + b)). (29) 

b is one of the three vectors that connect the nearest 
neighbors of the honeycomb lattice. The elements of the 
connectivity matrices T and K represent the uniform and 
Kekule hopping integrals between the nearest-neighbors, 
respectively. By rotating ri in the second term into t^ 
then, we find that the energy of the Dirac-Fermi sea in 
presence of a superconducting Kekule order parameter 
K is given by the sum of the energies of the two copies 
of the Dirac-Fermi seas for the spinless fermions: one 
in presence of the Kekule hopping pattern +K^ and the 
other in the pattern —K. We have therefore computed 
the energy /(a) for the single copy as a function of the 
parameter a at various values of the amplitude |A|. The 
typical result is depicted in Fig. 2. The function /(a) 
may be shown in general to be even, and periodic with 
the period 27r/3, which refiects the rotational symmetry 
of the honeycomb lattice. The computation shows that 
its absolute minimum is always at a = 0, in agreement 
with the recent work |20| . and the maximum at /(7r/3) — 
/(tt). We then find that 2/(7r/2) < /(O) + /(tt), as long 
as |A| < 2.725. The transition from the semimetallic 
phase is therefore into the superconducting Kekule phase 
with a = tt/2, which we therefore name p-Kekule. For 
|A| > 2.725, deep within the superconducting phase, we 
find a = solution to eventually become energetically 
favorable, with a discontinuous transition between the 
s-Kekule and p-Kekule superconductors. 



FIG. 2: The energy per site f{a) of the Dirac-Fermi sea of 
spinless fermions, hopping between nearest neighbors of the 
honeycomb lattice with the Kekule hopping amplitude 1 -|- 
jA|cos(Q ■ {x + y) + a) for |A| = 1, as a function of the 
parameter a. The precise type of the Kekule superconducting 
order depends on the sign of the combination /(O) + /{tt/S) — 
2/(7r/6). The points are the computed values, and the red line 
is our best fit -1.018 - 0.0175 cos(3a) -I- 0.000248 cos(6a) -I- 
O(10~^ cos(9q)). This implies the p-Kekule order (see the 
text). The transition into the s-Kekule state at |A| = 2.725 
essentially corresponds to the change in sign of the second 
harmonic of this function. 



VI. OTHER SUPERCONDUCTING STATES 

Let us also recognize the other gapped superconducting 
states, as the possible mass-terms that anticommute both 
with the Dirac Hamiltonian Hjj and with the number 
operator N: 

a) the standard s-wave superconductor with the on-site 
pairing, 

(*^ [(ri cos 4> + T2 sin (j)) (Exjq (g) 17073]*) , (30) 

which is translationally invariant, even under the valley 
and/or sublattice exchange, but odd under the exchange 
of spin labels (spin singlet). 

b) the f-wave Q 

(*'I'[(t-i cos(I) + T2 sin0) (g) ct® 17075]*), (31) 

which is translationally invariant, even under the sub- 
lattice and spin exchanges (spin triplet), but odd under 
valley exchange. 

c) the p-Kekule state discussed in the previous section 
written explicitly is 



(*n( 



Ti cos ( 



T2 sm ( 



>a(g)i7i72]*). (32) 



One can further construct all the gapless (hidden) con- 
densates, as 



ii^HlHo'^), 



(33) 



where M is a mass-matrix for any of the above gapped 
superconducting states. Choosing the matrix M to cor- 
respond to the s-wave superconductor yields the origi- 
nal hidden order of ref. (4). Since the Hermitian ma- 
trix iMHd by construction then anticommutes with the 




Kekule - SC 



FIG. 3: The schematic T — phase diagram in the model 
with the on-site ([/) and the nearest-neighbor (V) attractions. 
At the boundary between the two superconducting states the 
order parameter acquires the larger 0(5) symmetry, as the 
energy becomes invariant under the rotations of the Kekule 
into the s-wave order with the relative phase of 7r/2. 



Dirac Hamiltonian Hd while at the same time being pro- 
portional to it, its addition to Hu will effectively only 
renormalize the velocity of the Bogoliubov excitations, 
as manifest in Eq. (27). 

Of course, one can imagine many other matrices that 
do not commute with the number operator, and will 
therefore represent some superconducting order, which 
nevertheless do not fall into any of the categories listed 
above. These fail to anticommute with the Dirac Hamil- 
tonian, and as such neither gap out, nor increase the 
velocity of the Dirac fermions. Development of these or- 
der parameters would not therefore be particularly ener- 
getically advantageous, which is the reason behind their 
omission here. 

One can expect the gapped superconducting states to 
compete in the phase diagram for attractive interactions 
in a close parallel with the competition between insula- 
tors when the interactions are repulsive [l|. As an il- 
lustration, in Fig. 3 we present the mean-field phase 
diagram in presence of both the on-site and the nearest- 
neighbor attractions at half filling. In analogy with the 
insulating case, there is a discontinuous transition be- 
tween the ordered phases, whereas the transitions out of 
the semimetallic phase may be expected to be continuous 
|13 . Il3| . One novel feature is that because of the [/(I) 
symmetry in the superconducting phase the matrices rep- 
resenting different superconducting states can be chosen 
so as to anticommute. Consider the above p-Kekule state 
with the phase (j) — and the spin axis along z-direction, 
for example. Choosing the uniform s-wave state with 
(j) ~ tt/2 makes the two representative mass-matrices an- 
ticommuting, so that right at the boundary between the 
two phases the system acquires a larger symmetry 0(5). 



Adding the third axis for the second-nearest neighbor 
attractive interaction introduces then a region of the f- 
wave order, with the discontinuous transitions between 
any two of the three phases. Interestingly, there is a 
unique anticommuting f-wave state, with (j) — tt/2 and 
with the spin axis along z-direction, that may be added 
to the above combination of the already anticommuting 
Kekule and s-wave order parameters. At the point in the 
phase diagram where the three phases would meet, an 
even larger, 0(6), symmetry emerges. 

If the preferred non-uniform superconducting state is 
the s-Kekule, on the other hand, the possible anticom- 
muting states are again the uniform s-wave and f-wave 
condensates, but this time both with the same phase as 
the one of the s-Kekule state. 



VII. TOPOLOGY AND DEFECTS 

In the ordered phase, the mass-matrix for the Kekule 
order parameter in Eq. (26) lives on the Si x 5*2 tar- 
get space, but with opposite points identified. In other 
words, the order parameter space is the product of S2 for 
the spin direction and half of Si for the superconducting 
phase, which is equivalent to 6*3, the sphere in four di- 
mensions. That this is indeed the target space becomes 
clear upon recalling that the minimum with n — 0, the 
mean-field free energy in Eq. (21) depends only on the 
mass m^ — | A p -I- 1 A^ | ^ , where the two complex order pa- 
rameters A and A-|- are constrained only by the condition 
that m is the solution of the gap equation. 

Since the first and the second homotopy groups of S3 
are trivial, 7ri(S'3) = 7r2(S'3) = 1, there are no stable 
topological defects, and the massless fluctuations in the 
ordered phase should be correctly described by the 0(4) 
non-linear sigma model [21|. We therefore do not ex- 
pect a true finite temperature phase transition from a 
semimetal into the Kekule superconductor, but only a 
crossover when the superconducting correlation length 
^ ex exp(cTOo/T), with c as a (non- universal) numerical 
constant, reaches the size of the sample ^2\. 

A reduction of the rotational symmetry would change 
the target space and allow stable vortex excitations. Let 
us consider the case of a possible easy plane first. Such an 
anisotropy may be introduced most simply by placing the 
Kekule superconductor into a magnetic field. The Zee- 
man term representing the coupling of the magnetic field 
to the electron spin will be proportional to the genera- 
tor of rotations along the direction of the magnetic field. 
To <Si cr^ <^ I, for example. Since this matrix commutes 
with the Dirac Hamiltonian and with the third term in 
the Kekule mass-matrix M in Eq. (14) that is propor- 
tional to A, while it anticommutes with the two other 
terms in M that are proportional to A^., the minimiza- 
tion of the energy in the presence of Zeeman coupling is 
formally equivalent to the problem of Neel ordering in 
graphene catalyzed by the magnetic field [l^, 1231 ■ The 
result is that A = 0, since that way the Kekule mass 



matrix anticomniutes with the Zeeman term. The way 
to understand this physicahy is to reahze that in such 
a state the spins of paired electrons are all orthogonal 
to the magnetic field, so it becomes easier for them to 
tilt and provide a finite magnetization in the field direc- 
tion. The minimum condition n = then translates into 



|A-|-| = |A^|, the same as without the magnetic field, but 
without a further constraint on the phases of the two 
complex order parameters. The s-Kekule mass-matrix in 
Eq. (14) for such an easy plane may be then rewritten 
differently as 



J 



M 



TTlQ [T2 cos ■ 



Ti sm ■ 



-) <8) (o-2COS 



H- h 



(71 sm ■ 



-)«)7o- 



(34) 



r 



The target space with the easy plane anisotropy is thus 
S*! X 5*1 , with the factors corresponding to the two phases 
^"t- and 04,. Since the first homotopy group of Si is non- 
trivial, 7ri(S'i) = Z, there are different types of topolog- 
ically distinct vortex excitations. For example, winding 
just one of the phases by 2n causes both {(j)f + (j)^)/2 and 
(</)-[- — 04,) /2 to change from zero to tt, i. e. both the first, 
phase term, and the second, spin-axis term in the above 
matrix make half a circle. This is sometimes referred to 
as "half- vortex" [2J|. On the other hand, winding both 
the phases 0-|- and (j)^ in the same sense by 27r leaves the 
angle of the spin-axis intact, and produces the standard 
full vortex in the superconducting phase. Finally, wind- 
ing the two phases (1)-^ and 0^ in the opposite sense by 27r 
produces a third type of vortex, this time in the direction 
of the spin-axis only. 

An easy axis, on the other hand, is introduced by 
the spin-orbit coupling, for example. Consider adding 
a weak perturbation to the Dirac Hamiltonian propor- 
tional to T3 (g) CTa (g) J7172, which in our representation 
corresponds to the third component of the spin-triplet 
version of the time-reversal symmetry breaking mass, in- 
troduced by Haldane [g^ and discussed in the context 
of spin-orbit interaction in graphene by Kane and Mele 
[271 . The presence of such a term would again select the 
piece of the Kekule mass matrix that anticommutes with 
it, but due to the t^ matrix in the first, Nambu's factor, 
this now implies that A-j- = A^ = 0. The s-Kekule matrix 
assumes the form as in Eq. (26), with 9 = 0. The target 
space therefore in this case reduces to the usual Si , with 
only the standard vortices as the topological excitations. 
The internal structure of such a vortex has been studied 
in ref. [3]. 

Finally, it should be understood that we discussed 
the breaking of the rotational symmetry in the s-Kekule 
state for simplicity only, and that everything said applies 
equally to the p-Kekule state as well. 



VIII. SUMMARY AND DISCUSSION 

To summarize, we introduced the non-uniform super- 
conducting state on graphene's honeycomb lattice, and 
argued that it is the mean-field solution of the simple 



model with nearest-neighbor attraction. The order pa- 
rameter for this state lives on the bonds of the lattice, 
and forms the Kekule lattice with the period three. Com- 
petition between different such Kekule superconducting 
states, as well as between the Kekule and the other pos- 
sible gapped and gapless superconducting states was dis- 
cussed. The Kekule superconductor has a spin-triplet 
order parameter, which lives on the surface of the S3 
sphere. Target spaces for the order parameter and the 
topological defects in presence of some simple symmetry 
breaking terms were determined. 

A Kekule insulator which breaks the translational in- 
variance of the honeycomb lattice has been previously 
proposed and discussed in literature [lallal- However, 
it appears that this state is not the ground state of the 
simplest model with only the nearest-neighbor repulsion, 
since there is an energetically superior charge-density- 
wave that breaks the sublattice exchange symmetry avail- 
able. This should be contrasted with the situation for 
the attractive interactions considered here, where the 
competing superconducting state is the gapless super- 
conductor, which we argued should have a higher en- 
ergy. It was argued recently, however, that the Kekule 
insulator does become the mean-field ground state when 
there is a balance between the nearest-neighbor and the 
second-nearest-neighbor components of the repulsive in- 
teractions [2g |. 

We described here only the problem at half-filling in 
detail, where a finite interaction is needed to cause the su- 
perconducting transition. At a finite chemical potential 
the density of states at the Fermi level also becomes fi- 
nite, and there is the usual BCS instability at an infinites- 
imal attraction. For small deviations from the half-filling, 
however, the symmetry of the superconducting state is 
essentially determined by the solution at the Dirac point. 
For the nearest-neighbor attraction as the dominant com- 
ponent of the interaction one should therefore expect the 
non-uniform Kekule state we discussed to persist at a fi- 
nite doping as well. As long as the Fermi surface around 
the Dirac points stays circular [2g| the states with the 
momenta Q + q and Q — q may both be near the Fermi 
surface and be paired up by the Kekule order parame- 
ter with the momentum 2(5, essentially the same way as 
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right at half-filling. 
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X. APPENDIX 

Here we determine the the susceptibilities for the hid- 
den and Kekule orders used at the end of sec. IV, evalu- 
ated over the whole Brillouin zone. For the hidden order 
the energy per site may be written as 



^-iv-^T.\fCm\Af+oi\Ar) 



'V 2N 



where 



m= E 



ik-hi 



(35) 



(36) 



1=1,2,3 



-l-El/(fc)l -0.786 



2N 



(37) 



in agreement with [20| ■ This yields the value of V^ cited 
in the text. 

For the critical interaction for Kekule order we need 
the energy as a function of the Kekule mass m to the 
leading order. Diagonalizing the six-dimensional matrix 
given by Weeks and Franz 20] and summing over the 
reduced Brillouin zone for the Kekule lattice we find 



^(to) 

2N 



(A 

^2V 



0.727)m2 + O(|mp). 



(38) 



Note that the the electronic susceptibilities for the hidden 
and Kekule orders are rather close numerically, and the 
Kekule state wins mainly due to the geometrical factor 
of two in the first term in the mean-field energy. 
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